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Abstract 

We present measurements of the analyzing powers A y and iT\\ for proton-deuteron scattering at 
E c .m. = 432 keV. Calculations using a realistic nucleon-nucleon potential (Argonne V18) are found 
to underpredict both analyzing powers by 40%. The inclusion of the Urbana three-nucleon 
interaction does not significantly modify the calculated analyzing powers. Due to its short range, 
it is difficult for this three-nucleon interaction to affect A y and iT\\ at this low energy. The origin 
of the discrepancy remains an open question. 
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Few-body systems provide a fundamental testing ground for nuclear interactions. Com- 
parisons of measured three-nucleon scattering observables to theoretical calculations allow 
stringent tests of the underlying nucleon-nucleon (NN) and models of the three-nucleon (3N) 
interactions. Past studies have found that rigorous calculations utilizing realistic NN po- 
tentials underpredict by 25-30% the measured analyzing power A y in n-d scattering at low 
energies - a surprising discrepancy which has been dubbed the ll A y (9) puzzle" (see Ref. 
and references therein). Precise calculations for p-d scattering []2]|| below the deuteron 
breakup threshold including the Coulomb potential rigorously have recently become possi- 
ble to perform. These calculations find that a similar underprediction exists here for 
A y and also iTu. The agreement for other observables, including cross sections, T 20 , T 2 \, 
T 2 2, and n-d scattering lengths is generally excellent. The analyzing powers A y and iTu 
in N-d scattering are known to be very sensitive to the NN potential in the 3 Pj waves. It 
has been suggested that the potential in these waves may not be known to the necessary 
precision at low energies [§,0. The possibility that the underprediction is due to 3N force 
effects has also been considered This paper investigates the roles of NN and 3N force 

effects on A y and iTu for N-d scattering; we do not consider here other possibly important 
effects such as relativistic corrections or subnucleonic degrees of freedom. We do note, how- 
ever, that the Mott-Schwinger interaction, a long-ranged electromagnetic effect, has been 
recently shown not to be responsible for the discrepancies in A y and iTu H- 

Measurements of both A y and iTu are useful, as these observables depend on different 
combinations of phase shift and mixing parameters (principally in P- waves), and they also 
have been shown to be sensitive to different combinations of the 3 Pj NN interactions ||. 
The majority of the data on these analyzing powers has been obtained in the vicinity of the 
deuteron breakup threshold: n-d A y data exist for E cm > 2 MeV PJlOl p-d A y data exist 



for E c . m . > 0.53 MeV [TT|-|T3[, and p-d iT n data exist for E cm _ > 1.7 MeV |T2|,|T3[ (there 
are no n-d iTu data). It is desirable to determine these observables at lower energies, as 
the influence of higher partial waves is strongly reduced, and the dominant S- and P-waves 
can be investigated with more confidence. In addition, as we will show below, P-wave N- 



d scattering at low energies is almost entirely determined by the asymptotic part of the 
three-nucleon scattering wave function and the NN interaction. Under these conditions 
the connection between the measured observables and the underlying interactions is greatly 
simplified. The calculated analyzing powers A y and iTn are mainly determined by the j- 
splitting of the P-wave N-d phase shifts and the Ey 2 and e7, 2 mixing parameters. Due to the 
angular momentum barrier at low energies, these observables are very small and difficult to 
determine experimentally. This paper reports measurements of A y and %T\\ for p-d elastic 
scattering at E c , m , = 432 keV. These data are at the same energy as our previously-reported 
T 2 o and T 2 i measurements 0. 

These experimental results are compared to calculations utilizing the Pair- Correlated 
Hyperspherical Harmonic (PHH) basis |14| to construct the scattering wave function, and 



the Kohn variational principle to determine the scattering matrix elements . In addition, 
we present calculations using an "optimized" Born approximation |4| for the peripheral 
partial waves. The calculations have been done using the AV18 potential ]I5[ and with 
AV18 plus the 3N interaction of Urbana (UR) [TjJ. It has been shown in Ref. [|IJ that other 
high-quality NN potentials such as Bonn or Nijmegen predict essentially the same n-d A y 
just below the deuteron breakup threshold, so we would not expect our conclusion to change 
if these potentials were used. 

The measurements were performed using polarized proton and deuteron beams from the 
atomic beam polarized ion source |17j at the Triangle Universities Nuclear Laboratory. The 
deuteron beams were accelerated to E^ = 1.3 MeV using the FN tandem accelerator, and 
then directed into a 62-cm diameter scattering chamber. Proton beams were accelerated 
into a 107-cm diameter scattering chamber at E p = 650 keV using the minitandem accel- 
erator JIB] and chamber bias voltage ||19|| . Thin hydrogenated or deuterated carbon targets 
were utilized which consisted of approximately 1 x 10 18 and 1.5 x 10 18 hydrogen isotope 
and carbon atoms/cm 2 , respectively. The beams lose ~ 10 keV in these targets, with an 
average energy of E c m = 432 ± 1 keV for both measurements. The use of thin targets is 
very important at low energies for minimizing energy loss and straggling effects. 



The proton beam polarization was determined using the 6 Li(jo, 3 He) 4 He reaction in a 
polarimeter |2(J located at the rear of the scattering chamber. The polarization was mea- 



sured several times throughout the measurements at an incident proton energy of 450 keV 
by lowering the chamber bias voltage. The proton polarization was found to be constant 
within ±3% throughout the measurements; the systematic error in the proton polarization 
is estimated to be ±4%. Deuteron beam vector polarization was determined online via the 
12 C(d,p) reaction in a polarimeter located behind the scattering chamber. The effective iT\\ 
for this reaction at Ed = 1.3 MeV has been calibrated relative to the 3 He(d,p) reaction 



in another polarimeter at Ed = 12 MeV [21 1. The absolute uncertainty in the deuteron 



beam polarization is estimated to ±3%. For both beams the data were taken with the spin- 
quantization axis perpendicular to the reaction plane, using two spin states with p z ~ ±0.7 
for the proton beam; and pz ~ ±0.55, pzz ~ for the deuteron beam. The spin states were 
cycled approximately once every second, in order to minimize the effects of slow changes in 
beam position, target thickness, or amplifier gain. 

Scattered deuterons and protons were detected in coincidence using two pairs of silicon 
surface barrier detectors placed at symmetric angles on either side of the incident beam. The 
angles of the detectors were set to observe either protons or deuterons in the more forward 
detectors in coincidence with deuterons or protons detected in the more backward detector 
on the opposite side of the beam. Histograms of the time difference between the fast timing 
signals from each coincident pair of detectors were stored for each spin state. Dead-time 
corrections (< 3%) were determined by sending test pulses to the detector preamplifiers 
with time delays adjusted to give distinct peaks in the time spectra. The time resolution 
for the coincident proton-deuteron peaks was ~ 10 ns, with backgrounds < 3%. The ana- 
lyzing powers were determined from the counts in the coincident peaks, after correction for 
background, dead time, and the number of incident particles (determined by beam-current 
integration). It should be noted that the coincidence technique is essential for measuring 
the small analyzing powers A y and %T X \ with these targets, as the elimination of carbon 
elastic-scattering events by the fast coincidence requirement allows proton-deuteron scatter- 



ing events to be counted at the high rate required to achieve reasonable statistical accuracy. 
The results for A y and iTu are shown in Fig. [1]. The error bars include contributions from 
statistics and background subtraction, but not the absolute beam polarization. 

The theoretical method has been described previously |||§; it can be applied equally 
well to n-d as well as p-d scattering, and realistic NN and 3N potentials can be used without 
difficulty In the present calculations, scattering waves with orbital angular momentum up 
to L = 4 have been taken into consideration. At this energy, the differential cross section, 
A y , and iTu are determined almost entirely by waves with L < 1, while for T20, T21, and 
T22 L = 2 waves are also important. In particular, A y and iTu change by < 10~ 4 when 
phases with L > 1 are considered. In Fig. [I] the data are compared to the calculations using 
the AV18 potential and the AV18+UR potential. The corresponding P-wave and 4 Si/ 2 
phase-shift parameters are given in Table |I[ 

It is seen that both calculations underpredict the data by » 40%. The change in the 
calculated A y and iTu resulting from the inclusion of the 3N interaction is too small by an 
order of magnitude to explain the discrepancy. 

Relatively small changes in the N-d phase shift parameters can have large effects on 
the corresponding analyzing powers. In Ref. || it was found that the discrepancies in A y 
and iTu for E c _ m , = 1.67 and 2 MeV could be corrected by reducing the A P\/2 phase shift 
by 3.4% and increasing the absolute value of Ey 2 mixing parameter by 12%. Using the 
AV18+UR results for the other phase-shift parameters, we find that agreement with our A y 
and iTu results is optimized if the 4 Pi/2 and e^ 2 parameters given in Table | (column 2, 
in parenthesis) are replaced by 5.22° and —1.02°, i.e., the absolute values are reduced by 
1.6% and increased by 15%, respectively. The results for A y and iTu using these parameters 
are shown by the long-dashed curve in Fig. [I| It is important to note that these parameter 
changes affect the cross section by < 0.15%, and the other analyzing powers by < 0.0012. 
In particular, the good agreement observed previously with our T20 and T21 data at this 
energy |jj is not disturbed. While similar to the changes required at higher energies [|J, 
there are significant differences in the fractional changes required. We should point out 
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however that there is no reason to expect the percentage change required to be the same for 
different energies. 

We have also performed calculations using an "optimized" Born approximation, i.e., the 
procedure in which the second-order 7?.-matrix is estimated using the asymptotic part of 
the three-nucleon wave function as described in Ref. ||. The results for the 4 S^/2 and P- 
wave phase-shift parameters are given in Table | and compared to those obtained when the 
complete wave function is considered (full solution). It is seen that the Born approximation 
results are close to the full calculation for these partial waves. In the case of the A Sz/2 phase 
the Pauli principle prevents the three particles from being close to each other, while for the 
P-waves the centrifugal barrier is sufficiently high at these energies. These findings indicate 
that these partial waves are almost entirely determined by the asymptotic structure of the 
system. We also show in Table p] the results for the full solution and Born approximation at 
E c .m. — 2 MeV, where it is seen that the accuracy of the Born approximation for low partial 
waves is reduced. On the other hand we observe that the influence of the 3N force is small 
and of the same magnitude at both energies. 

The small effect on A y and iT n from including the UR 3N interaction is now clear. 
This interaction, which is based on two-pion exchange and includes a phenomenological 
repulsive short range term, requires the three nucleons to be close together in order to 
produce a significant effect. The likelihood for this situation is diminished by the diffuse 
structure of the deuteron which results from the small binding energy. For low energy p- 
d scattering in P-waves (or higher L values), the probability of finding three nucleons in 
close proximity is further reduced by the centrifugal and Coulomb barriers. We thus draw 
the important conclusion that 3N interactions based on two-pion exchange cannot produce 
significant changes in A y and iT n at low energies. 



Other choices for the 3N potentials, such as the Tucson-Melbourne |22|] or the Brazil |23 
models, give quite similar conclusions. Inclusion of other processes, such as n — p or p — p 
exchanges, involving heavier mesons and therefore shorter ranges, are expected to give still 
smaller corrections [0J. These findings thus indicate that new types of 3N interactions should 



be considered. One possibility is the inclusion of a spin-orbit 3N force which could signifi- 
cantly affect the N-d P- waves []2|| . One cannot exclude also the possibility that inadequacies 
in the NN interaction are responsible for the discrepancy. 

In summary, our measurements of A y and iT\\ at B m = 432 keV are significantly 
underpredicted by calculations utilizing the AV18 NN interaction. The inclusion of the UR 
3N interaction does not significantly change the theoretical calculations. We have shown 
it is difficult to identify a 3N interaction which could significantly change these analyzing 
powers at low energies, as they are mainly determined by long-ranged interactions. It would 
be of great interest to extend these comparisons to p- 3 He and n- 3 H scattering, where the 
3N force effects are expected to be larger, as the likelihood of finding three nucleons close 
together is enhanced by the tighter binding of 3 H and 3 He. 
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FIGURES 

FIG. 1. Experimental A y and iT\\ for p-d scattering at E c . m . = 432 keV (circles), along with 
theoretical calculations using the AV18 (solid line) and the AV18+UR potentials (short-dashed 
line). The long-dashed line results from modifying the 4 Pi/2 and Ey 2 phases as described in the 
text. 
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TABLE I. Some of the theoretical phase shifts and mixing parameters (in degrees) calculated 
at E c . m . = 0.432 and 2 MeV for the AV18 potential; the values in parentheses correspond to the 
AV18+UR potential. For both energies, the results for the "optimized" Born approximation and 
the full solution are reported in the columns labeled "Born" and "Full" , respectively. 



12 



0.0080 




50 100 

e c.m. (degrees) 



